pss header will be provided by the publisher 



Thermoelectric effects in STM tunneling through a monoatomic 
chain 

M. Krawiec* and M. Jatochowski 

Institute of Physics and Nanotechnology Center, M. Curie-Sklodowska University, pi. M. Curie- 
Sklodowskiej 1, 20-031 Lublin, Poland 

Received 9 September 2006 

Key words STM; atomic wire; tunneling 
PACS 68.37.Ef, 81.07.Vb, 73.40.Gk 

We study thermoelectric properties of the system composed of a monoatomic chain on a surface and ad- 
ditional electrode coupled to the chain, which can be an STM tip. In particular, we are interested in ther- 
mopower, electric and thermal conductance, Wiedemann-Franz relation and thermoelectric figure of merit, 
which is a direct measure of the usefulness of the system for applications. We discuss the modifications 
of the STM wire topography due to temperature gradient between the electrodes. Finally, we also make 
connection to STM experiment, in which the thermopower has been directly measured, showing different 
structure, not visible in topography spectra. 
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1 Introduction 

Since the invention of scanning tunneling microscopy (STM) [1] it became possible to tailor and ana- 
lyze small nanostructures on various conducting surfaces JU [3J. As the STM is a real space technique 
and is very sensitive to local atomic and electronic structures, it has been widely applied to study vari- 
ous surface reconstructions [3] and low dimensional structures, like single atoms 12E), islands ||4][5) or 
one-dimensional monoatomic chains [6|-[9|. While the STM is used to characterize various structures 
measuring topography or current-voltage characteristics, it is also suitable to get information on thermal 
properties of the system. For example, local thermopower S(r) of various surfaces has been measured with 
STM ifTOlfTTl . Thus the topography of the surface can be supplemented by its 'thermal' image. It turns out 
that the thermal images show features that are not accessible in the conventional STM topographic images. 
In particular, the maxima of the thermopower in the thermal image of the surface do not overlap with the 
tunneling current maxima IfTOlfTTl . The local thermopower of the surfaces observed in STM experiments 
IfTOlfTTl has been also studied theoretically [12] within the Tersoff-Hamann approach [13]. It was found 
that the magnitude of the thermopower depends on the logarithmic derivative of the local surface density 
of states at the Fermi level and does not show exponential dependence on tip-surface distance, in contrast 
to the tunneling current. 

The above examples show that thermoelectric properties are the source of information complementary 
to that obtained from other transport characteristics. In bulk systems, transport in the presence of electrical 
and thermal gradients is a well studied phenomenon. Recently, it became possible to study those effects, 
both experimentally and theoretically, in atomic-scale systems, like quantum point contacts lfl4l[T5l . quan- 
tum dots [16|-[21|, quantum dot superlattices [22], quantum wires with end atoms coupled to external 
leads [23 1, carbon nanotubes [24 25 1 or correlated multilayered nanostructures | 
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It is the purpose of the present work to study thermoelectric properties of quantum wire in STM geome- 
try. We shall concentrate on the electric and thermal conductance, thermopower and related quantities, like 
thermoelectric figure of merit which is a direct measure of the usefulness of the system for applications and 
Wiedemann-Franz ratio which signals breakdown of the Fermi liquid state. We also discuss the modifica- 
tions of the STM wire topography due to temperature gradient between the electrodes. To explore electric 
and thermal properties of the system we propose a model of tunneling between STM tip and surface and 
use equations of motion technique for Green functions. Rest of the paper is organized as follows: in Sec. [2] 
we introduce our model based on tight binding model and discuss some aspects of our procedure, and the 
results of our calculations are presented in Sec. [3] We end up with summary and conclusions. 



2 The model 

Our model system is described by the Hamiltonian 

Ae{t,s}kcr a in ija 

+ E ( v t c t + kCT c °- + H - c ) + E + h.c.) + e {u 4<,*° + h.c.) , (i) 

kcr jker ia 

and consists of N-atom metallic wire with atomic energies e w and hopping parameter t w between neigh- 
boring wire atoms. The wire is connected to the surface via parameters V; s , which we treat as a reservoir 
for electrons with energies e s k- STM tip is modeled by single atom with energy Eq attached to another 
reservoir (with electron energies e t k) via coupling V t . Tunneling between STM tip and the wire atoms is 
described by tunneling matrix element t^o- As usually, c\ (c\) stands for creation (annihilation) electron 
operator in STM lead (A = t), tip atom (A = 0), ith wire atom (A = i) and surface (A = s). Schematic 
view of our model system is shown in Fig. 1 . 




Fig. 1 Schematic view of our model system composed of N-atom wire on surface and STM tip. 



In order to calculate tunneling current J e and thermal flux Jq flowing from the STM electrode to the 
rest of the system we follow standard derivation ||2T1 and get 

J e = | J duT(uj) [f(ui - /i t ) - f(ui - n s )] , (2) 

J Q = \ J dujT(uj)(uj - m) [f{u) - m) - f{u) - /j, s )] , (3) 

where T(cu) = 2r* (uj) r \- (uj)Gj (w)Gq^ (w) is the transmittance of the system, J(lu) is the Fermi dis- 
tribution function, and fit (Ms) is the chemical potential of the tip (surface) electrode. (w) is the matrix 
element (connecting the tip atom with the jth wire atom) of the retarded (advanced) Green function G(oj), 
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the solution of the equation (luI — H)G(u>) = 1. The parameter r t (cj) = ^ k |Vt| 2 <5(o; — e t k) denotes 
strength of the coupling between the tip atom and STM electrode, while rf-(w) = |V s | 2 e lk ^ Ri_R ' j ' ) (5(a;- 
e t k) = T s (oj)sm(kpa(i — j))/(kpa(i — j)) is the coupling between wire atoms i and j via the surface. 
kp is the Fermi wave vector of the surface electrode, and a is a distance between neighboring wire atoms. 
In typical metals fcpa is of order of 4 — 5. 

For small bias voltages eV = /it — Ms - * and small temperature gradients ST = T t — T s — > one 
defines conductance G = — (e 2 /T)Ln, thermopower is given in the form S = — (l/eT)(Li2/Ln), and 
thermal conductance k = (1/T 2 )(L22 — L\ 2 /Ln). The linear response coefficients read 



/i J \ djji 



t 2 /■ , w 



i 22 = T J (kjT(cj)(u - Mt ) (^^^ ) . (6) 



3 Results and discussion 

Before the presentation of numerical results, we would like to comment on choice of the model parame- 
ters used in the present work. In numerical calculations we have assumed equal and energy independent 
coupling parameters (T^^u;) = r*^) and chosen r s = r* = T as an energy unit. The other parameters 
have been chosen in order to satisfy realistic situation in experiments. The hopping integral along the wire 
is t w = 2, the parameter connecting tip with underneath wire atom tjo = 0.1 and STM tip and wire atomic 
energies Eo = e w = 0. For example, taking T = 0.05 eV, we get t w = 0.1 eV and tjrj = 0.005 eV. Such a 
value of t,o gives tip-surface distance z;o = 6 A. Also temperature is measured in units of V (ks — 1) and 
for example, T = 1 corresponds to w 500 K. 

Figure 2 shows conductance G (top panels), thermal conductance k (middle panels) and thermopower 
S (bottom panels) of quantum wire. The temperature dependencies of those quantities for wires consisted 
of different number of atoms (indicated in the figure) are shown in left panels, while in right panels the 
same quantities are plotted as a function of number of atoms in the wire at different temperatures. In all 
the cases the STM tip is above first wire atom, and the tip-surface distance is equal to 6 A. 

The conductance G shows rather an expected behavior, i.e. at high T it goes like T _1 , while in the 
low temperature regime remains constant. Note, that it never reaches unitary limit (2e 2 /h) due to large tip 
surface distance. However there are differences in maximal values of G for different number of atoms N in 
the wire. To be more precise, the conductance shows even-odd oscillations with increasing N (see right top 
panel). Such a behavior is similar to that observed in the transport along wires l27l . Similar oscillations 
shows thermal conductance k, although they are less pronounced at higher temperatures (see right middle 
panel). However, it has different temperature dependence, it goes linearly with T at low temperatures, 
while at high T shows T~ 2 behavior. The thermopower S goes linearly to zero with T and behaves like 
T _1 in high temperature regime. It is almost always positive in the whole range temperatures, suggesting 
an electron nature of transport. Interestingly, it shows also even-odd oscillations with number of atoms in 
a wire but with opposite phase to G and k, i.e. while G and k show maxima for the wires composed of odd 
number of atoms, S has maximal values for even N. This can be easily understood, as the conductance G is 
sensitive to density of states (DOS) at the Fermi level Ep, while S reflects the curvature of DOS around it. 
In odd atom wire the DOS is large at Ep and has small curvature, while in even atom wire the situation is 
opposite. Note that the presence of the surface introduces asymmetry to the DOS around the Fermi energy, 
even if we assume wire energies to coincide with Ep Il28l l9ll. 
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Fig. 2 Conductance (top panels), thermal conductance (middle panels) and thermopower (bottom panels) as a function 
of temperature (left panels) and as a function of number of atoms in the wire (right panels). All the curves correspond 
to situation when STM tip is placed above first wire atom, and the tip surface distance is equal to 6 A. 



Topography of the wire consisted of N = 19 atoms at eV = 0.1 and J e = 10 5 (V = 0.5 mV and 
J e = 0.2 nA) is shown in Fig. 3. Dashed line corresponds to the situation when STM tip and the surface 




Fig. 3 Topography of the wire z(x) at eV = 0.1 and J e = 1CT 5 (V = 0.5 mV and J e = 0.2 nA) for the STM tip 
temperature T t =0.1 and different temperatures of the surface electrode T s (top panel). Corresponding thermopower 
S(x) is shown in the bottom panel. For comparison z(x) is also displayed, x and z are in units of the distance between 
neighboring wire atoms a. 



temperatures are equal (T t = T s = 0.1 (50 K)), while the solid one to the situation when surface is at 
elevated temperature T s = 0.5 (250 K). The reverse of the topography is clearly seen due to the temperature 
difference between STM tip and the surface. To be convinced that this is really the case, we have also 
shown thermopower S(x) (solid line in the bottom panel). Note that the maxima of S(x) do not coincide 
with the topography maxima calculated for Tt = T s (dashed line in the bottom panel) but they do with 
those obtained at different STM tip and surface temperatures (dashed line in the top panel). Thus indeed 
the reverse of the topography is caused by the temperature gradient between STM tip and surface. Similar 
differences between the topography and thermal images of the surface have been observed experimentally 
with STM HI- 
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The Wiedemann-Franz (WF) relation k/TG = 7r 2 /3e 2 describes transport in Fermi liquid bulk metals 
and in general is not obeyed in mesoscopic systems. In our system there are no Coulomb interactions and 
this relation is fulfilled at low temperatures, thus indicating Fermi liquid ground state. This can be also 
understood from temperature behavior of electric and thermal conductances, as k decreases linearly with 
T, while G remains constant. At high temperatures the WF relation is violated due to large fluctuations of 
G, which are of order of 2e 2 /h ||29l . 

Finally, we would like to comment on usefulness of such systems for potential applications, like thermo- 
electric power generators or cooling systems l22l . A direct measure of it is thermoelectric figure of merit 
Z = S 2 G/k. For simple systems it is inversely proportional to operation temperature, thus conveniently 
consider ZT, which numerical value is a measure of the system performance. In our system the value of 
ZT is always smaller than 1, indicating limited practical applicability. Also the thermopower, which is of 
order of 10~ 2 mV/K, leads to similar conclusions. 

4 Conclusions 

In conclusion we have studied thermal properties of monoatomic quantum wire in STM geometry, and 
show differences between topography and thermal images of the wire. In particular, the maxima of cal- 
culated thermopower along the wire do not coincide with those of topography spectra due to the fact that 
tunneling current is sensitive to the density of states, while thermopower is sensitive to the slope of the den- 
sity of states at the Fermi energy. We discussed also Wiedemann-Franz relation and found that is fulfilled 
at low temperatures indicating Fermi liquid ground state. 
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